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Abstract 

We design a simple ascending-price algorithm to compute a (1 -I- £)-approximate equilibrium in 
Arrow-Debreu exchange markets with weak gross substitute (WGS) property, which runs in time 
polynomial in market parameters and logl/e. This is the first polynomial-time algorithm for most of 
the known tractable classes of Arrow-Debreu markets, which is easy to implement and avoids heavy 
machinery such as the ellipsoid method. In addition, our algorithm can be applied in unknown market 
setting without exact knowledge about the number of agents, their individual utilities and endowments. 
Instead, our algorithm only relies on queries to a global demand oracle by posting prices and receiving 
aggregate demand for goods as feedback. When demands are real-valued functions of prices, the oracles 
can only return values of bounded precision based on real utility functions. Due to this more realistic 
assumption, precision and representation of prices and demands become a major technical challenge, 
and we develop new tools and insights that may be of independent interest. 

Furthermore, our approach also gives the first polynomial-time algorithm to compute an exact 
equilibrium for markets with spending constraint utilities, a piecewise linear concave generalization of 
linear utilities. This resolves an open problem posed in [^ . 
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1 Introduction 


The concept of market equilibrium is central in economics and captures fair, stable, and efficient outcomes 
in competitive allocation scenarios. The most prominent model to study market equilibria are exchange 
markets [1], which consist of a set of divisible goods and a set of agents. Each agent has an initial 
endowment of goods and a utility (preference) function over bundles of goods. Given prices of goods, 
each agent buys a most preferred bundle (called demand) that is affordable from the earned money. At 
equilibrium, market clears, i.e., demand meets supply. 

The computation of market equilibrium is a fundamental problem in economics and computer sci¬ 
ence [21ESI123135] • The challenge is to provide algorithms to compute an (approximate) market equi¬ 
librium efficiently. After more than a decade of research on this issue in theoretical computer sci¬ 
ence [71[51 fT3HT3l[2ni[^[^[^[32ll36lH5] , there is a fairly good understanding of tractable and intractable 
domains (assuming P ^ PPAD). For exchange markets, the tractable case is essentiall}0 given by markets 
with a weak gross-substitutes (WGS) property, where any increase in prices of a set of goods does not 
strictly decrease the demand of untouched goods. This includes markets with many popular and inter¬ 
esting classes of utility functions [33|, for example utilities with constant elasticity of substitution (GES) 
with 0 < p < 1, Cobb-Douglas utilities, linear utilities, or utilities with spending constraint^ [22ll42j . 
The only general approach to compute an approximate market equilibrium for this class of markets uses 
the ellipsoid method and convex feasibility programs [14j . 

Algorithms based on the ellipsoid method, although being polynomial time in principle, have a number 
of undesirable drawbacks. First, the ellipsoid method often tends to be slow in practice. Moreover, such 
algorithms often are not informative about the problem structure since they invoke a powerful algorithmic 
tool in a “black-box” fashion. In many algorithmic domains it has become an important line of research 
to avoid black-box application of the ellipsoid method and derive simpler, “combinatorial” algorithms 
that reveal inherent structure. Thereby, such algorithms also avoid possible overhead that comes with 
the application of broad algorithmic hammers and tend to be much faster. Despite extensive research [la- 
the search for a simple and direct polynomial-time algorithm for computing market equilibria has 
so far been unsuccessful. The only exceptions are a classic strongly polynomial-time algorithm for the 
case of Cobb-Douglas utilities [26], and recent combinatorial algorithms for linear utilities [23[I24| . 

A prerequisite for all these efficient algorithms is the entire description of the market including number 
of agents, their utility functions and initial endowments. Obtaining this description is a highly non-trivial 
task; an entire theory of revealed preferences [371I401HT] was developed to study how to infer market 
parameters from observed prices and buying patterns. Further, sometimes agents’ preferences can be 
complicated, and their utility function may not have any succinct representation. This gives rise to the 
following question: Can we design efficient algorithms that are oblivious to market parameters? In other 
words, are there efficient algorithms for unknown markets? 

Even beyond the economic interpretation, algorithms for computing market equilibria have wide ap¬ 
plicability. For example, proportionally-fair allocations, which result from market equilibria, are widely 
used in the design of computer networks [33]. Recent applications also include energy-efhcient schedul¬ 
ing [3011311 and fair allocation of indivisible resources [17[I18]. In the latter application, a combinatorial 
polynomial-time algorithm developed for linear markets [33] is extended to a problem, for which no con¬ 
vex programming formulation is known. This further highlights the importance of designing algorithms 
without the ellipsoid method that are simple and flexible to adjust to applications in related domains. 

In this paper, we design simple ascending price algorithms for WGS exchange markets. Under the 

^The other tractable cases are CES utilities with — 1 < p < 0 m, and piecewise linear concave utilities provided that 
the number of goods are constant m- 

^Spending constraint utilities are a piecewise linear concave generalization of linear utilities, which satisfy many desirable 
properties and have many additional applications. 
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standard assumptions that the most preferred bundles are unique and continuous in all prices, our algo¬ 
rithm computes an (1-|-e)-approximate market equilibrium in polynomial time, i.e., in time polynomial in 
market parameters and log(l/e). In markets with linear and spending constraint utilities the uniqueness 
and continuity assumptions do not hold. Thus, even with tie-breaking rules that imply the WGS prop¬ 
erty, the existing algorithms for general WGS markets are not applicable. Instead, our algorithm can be 
adapted when using an appropriate tie-breaking rule for demands. We improve over the FPTAS [22j in 
terms of running time from 1/e^ to log(l/e). More importantly, the log(l/e)-dependence of the running 
time allows us to convert the approximate equilibrium to an exact one in polynomial time when given 
all utility parameters. We obtain the first polynomial-time algorithm for computing an exact market 
equilibrium for spending constraint utilities and settle an open question raised by |24j . 

All our results are achieved by a unifying framework that works directly on the price vector. It uses 
simple binary search to identify suitable multiplicative price updates for subsets of goods. As such, our 
approach is extremely easy to implement and, in particular, avoids black-box use of the ellipsoid method. 
Furthermore, by working on the full generality of exchange markets, our algorithms can be applied to the 
simpler case of Fisher markets [6]p| . 

Moreover, to compute approximate equilibria our algorithms do not require access to an explicit 
description of the utility functions and endowments of the agents. Instead, the algorithms even work in 
what we term unknown markets - the number of agents, the agents’ endowments and utilities all remain 
unknown. As in the case of tatonnment algorithms [3l-[5lllf) t fT H ll6 p 39] . we assume these parameters can 
only be queried implicitly via aggregate demand queries. In such a query, we present a non-negative 
price for each good. For each agent, the oracle translates the endowment of each agent into money and 
determines a utility-maximizing bundle of goods for the money. It then aggregates the demands for 
each good and returns the vector of aggregate demands as answer to the query. Note that for many 
WGS markets these demand oracles can be easily implemented, e.g., for CES utilities there are even 
closed-form formulas. For suitable tie-breaking, we can even apply the algorithms to unknown markets 
with linear and spending constraint utilities, which have non-continuous demand. Perhaps surprisingly, 
our tatonnement-style algorithms succeed to detect implicitly all relevant information and changes in 
preferences in an unknown market setting. This requires to overcome a number of technical challenges 
and introduction of new techniques that we discuss in more detail below. 

1.1 Model and Notation 

Exact and Approximate Market Equilibrium In an exchange market there is a set A of n agents 
and a set G of m goods. Each agent i has an initial endowment Wij S ]R>o of good j. We denote the 
total endowment of good j by Wj = J2i£A'^ij assume w.l.o.g. that wj = l,Vj G G. An allocation 
X = {xij)i£A,jeG is an assignment of goods to agents such that Xij > 0,Vi G A,j G G, and = 1- 

Each agent i has a utility function Uj(xj) which specifies the value agent receives from his bundle Xj. 
An (exact) market equilibrium is a pair (x, p), where x is an allocation and p = {pj)j^G is a vector of 
non-negative prices pj >0. In a market equilibrium, each agent obtains a budget of money by selling 
his endowment at the given prices. Then Xj represents a utility-maximizing bundle of goods that he can 
afford to buy at the current prices for his budget. We call such a bundle a demand bundle of agent i. In 
addition, we say Xj = J2i ^ij is the demand for good j in allocation x. 

Defiuitiou 1.1. A bundle Xj is a demand of agent i at prices p ifui{'Xi) = max{uj(yi) | y* G K>q, p'^yj < 
p^Wj} . A pair (x, p) is a market equilibrium if (1) :x.i is a demand of agent i at p, and (2) ^ij = J2i '^ij 

^Fisher markets are a subclass of exchange markets where buyers and sellers are different agents. Buyers bring money to 
buy goods while sellers bring goods to earn money. 
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for each good j £ G. If a pair (x, p) is not a market equilibrium, the excess demand Zj = Xj — 1 is non-zero 
for at least one good j £ G. 

Let us consider a concept of strong approximate market equilibrium M- It relaxes only the market 
clearing constraint but not the condition that Xj is a demand for each agent i. 

Definition 1.2. A pair (p,x) is a ;U-approximate equilibrium {fj, > 1) if (1) for each agent i, Xj is a 
demand of agent i at prices p, and (2) ^ fof ^CLch good j. 

In Fisher markets [6], agents are divided into buyers and sellers. Each buyer i comes with an initial 
endowment of money Bi. Each seller i has an initial endowment of goods Wj. Each buyer has no value 
for money and is only interested in buying goods, and each seller is only interested in obtaining money. 
Eisher markets are a special case of exchange markets when we interpret money as a separate good. 

Utility Functions Algorithms for computing market equilibria rely on structural properties of the 
utility functions. A natural class are linear utilities when each agent i has non-negative values Uij > 0 
for each good j £ G, and ?ri(x) = ^ generalization, we consider spending constraint 

utilities [22], where the utility derived by agent i from good j is given by a piecewise linear concave (PLC) 
function fij. The overall utility of agent i is additively separable among goods, i.e., Wi(x) = J2j£G 
Each fij is a PLC function with a number of linear segments. Each segment k has two parameters: the 
rate of utility Uijk per unit of good derived on segment k and the maximum fraction Bijk of budget that 
can be spent on segment k. All Bij^ are strictly positive, and concavity implies Uijk > Here we 

assume all UijkS are integers, all Bijk,Wij's are rational numbers and the whole input can be represented 
in no more than L bits. Markets with spending constraint utilities have an equilibrium composed of 
rational numbers under mild sufficiency conditions (see Section [3] for details). 

More generally, we consider non-decreasing utility functions that generate markets with the weak 
gross-substitutes (WGS) property - when we increase a price, the demand for goods with untouched 
prices does not strictly decrease. For general WGS markets, we will assume that all demand bundles 
of agents are unique. Prominent examples are markets with Cobb-Douglas utilities Ui{:x.) = Y\x^2^ or 

constant-elasticity-of-substitution (CES) utilities rtj(x) = with 0 < p < 1. Note that 

even if all utility and endowment parameters are rationals of finite size, then demand bundles and the 
market equilibrium in such markets might involve irrational numbers. In this case, we are interested in 
approximate market equilibria, and our prices will use a prespecified precision depending on the desired 
approximation factor. 

Oracles Our algorithm queries demands for the agents by publishing prices p. Then an oracle returns 
the total demand Xj for each good f £ G. It assumes each agent can sell his initial endowment at the 
given prices and then requests a utility-maximizing bundle of goods for the money he has available. For 
ease of notation, given any price vector p, let 0(p) denote the surplus vector s = (si,..., Sm) for the 
return of the oracle, where Sj = pjZj is the surplus (in terms of money) of good j £ G. In other words, 
assuming we publish p, the oracle returns an excess demand vector z = (zi,..., Zm), ther0 0(p) = p • z. 

In general WGS markets the surplus vector might contain irrational values. Thus, we use an approx¬ 
imate demand oracle 0(p,//), which is a blackbox algorithm that takes any price vector p and positive 
rational p as input. It returns a surplus vector s such that |si — 0(p)| < p holds for every good i. Note 
that for many WGS markets the demand oracle can be implemented very quickly - for GES utilities there 
are even closed-form expressions for the demand of each agent for each good as a function of prices, utility 

"^We define the surplus vector as price times excess demand since it satisfies ^ - Si = 0. This invariant is useful in design 
and analysis of our algorithms below. 
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and endowment. More generally, we assume that the oracle can be implemented in time polynomial in 
the input size and log(l///). This standard assumption for demand oracles has also been used in previous 
work [TBIITI] . 

In linear and spending constraint markets, a major challenge for an algorithm in unknown markets is 
non-uniqueness of demands. Here an oracle needs to do tie-breaking between several different bundles of 
goods that yield the same maximum utility for an agent. Ideally, it should satisfy the following properties: 
(1) The output demand is always deterministic and unique. (2) If p are equilibrium prices, the output 
demand equals supply for every good. (3) The oracle can be implemented in time polynomial in the input 
size. Based on these criteria, we use a demand oracle that yields demands minimizing the £ 2 -iiorm of the 
surplus vector. More formally, for prices p our oracle returns a set of demand bundles for the agents such 
that J2i is minimized, where Si is the surplus of good i. Such a tie-breaking rule was introduced by m 
and it has been later used in several market algorithms, such as |23p 24 j H2]. This oracle satisfies all three 
properties. Furthermore, if utilities, endowments and prices are all given as integers with a number of 
bits polynomial in m and L, we can represent every surplus Si also with a number of bits polynomial in 
m and L. Hence by setting /r sufficiently small, we can convert an approximate demand oracle O into an 
exact oracle in polynomial time. Therefore we will assume that in spending constraint markets we are 
equipped with an exact demand oracle O instead of an approximate one. 

1.2 Results and Contribution 

We present simple ascending-price algorithms that converge to market equilibrium in WGS and spending 
constraint markets. Our algorithm for WGS markets in Section [2] converges to a (1 -|- e)-approximate 
equilibrium in time polynomial in m, log(l/e), and other market parameters. We present the first 
algorithm which is not based on the ellipsoid method for this general class of markets. Furthermore, 
the number of agents, the agents’ endowment, and utilities all remain unknown. We only query a global 
demand oracle that provides aggregate demands for goods at given prices. This information is then 
used to increase prices of a selected set of goods whose demand is more than their supply. In Section [5] 
we ignore for simplicity all precision and representation issues to highlight the general proof technique. 
The complete analysis with all details of our algorithm is provided in Appendix [Al where we specify 
in advance a precision for prices and rely on approximate demand oracles whose output is within our 
desired bit precision. One can view our algorithm as a form of tatonnement. It improves over previous 
approaches [13] by exponentially decreasing the dependence of the running time on 1/e and by working 
directly with the unknown market without transformations like adding auxiliary agents. 

Next, we apply our approach in Section [3] to spending constraint utilities - a piecewise linear concave 
generalization of linear utilities, which has many additional applications due to its natural diminishing 
returns property. Since these markets have demands that are non-continuous in the prices, we cannot 
directly apply our algorithm or other previous algorithms for WGS markets. Instead, we adjust our 
approach to implicitly capture the non-continuity events when using a global demand oracle with suitable 
tie-breaking. When all parameters are represented by at most L bits, our algorithm computes even an 
exact market equilibrium in time polynomial in m and L. All prices and demands occurring during 
the algorithm require a bit precision polynomial in m and L. It first computes a (1 -|- e)-approximate 
equilibrium using a precision that is polynomial in m, L and log(l/e). The exact demands returned by 
the demand oracle have the same precision. For a small enough e (using only polynomial bit length), we 
can then use a rounding procedure to turn it into a price vector of an exact market equilibrium. 

Note that our algorithm requires only access to a suitable demand oracle to compute an approximate 
equilibrium. However, in contrast to WGS markets, for spending constraint markets the oracle uses global 
tie-breaking. To obtain an exact equilibrium, our final rounding procedure relies on full information about 
the utilities. Thus, we obtain in polynomial time an approximate equilibrium in unknown markets (with 
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global tie-breaking) and an exact equilibrium with full information. This represents the first polynomial¬ 
time algorithm to compute an exact market equilibrium for spending constraint utilities and settles an 
open question raised by [21] • An important open problem is to construct efficient tatonnement algorithms 
that avoid global tie-breaking and full-information rounding. 

We use and extend ideas of algorithms for linear markets with full information I2I1IMI- These ideas 
were also used in spending constraint markets |22j to compute a (1 -|- e)-approximate equilibrium in time 
polynomial in m, L and 1/e. Roughly speaking, they are intimately tied to linear and spending constraint 
utilities, where they work on the agents’ side and increase prices until structural changes occur in the 
optimal bang-per-buck relations. The progress towards equilibrium is measured in the reduction of the 
^ 2 -norm of surplus. Our approach works for all WGS markets, in which in general do not exhibit such 
structural events that can be used for analysis. This turns out to be much more demanding, and we 
develop an approach that works entirely on the goods’ side, based only on prices and aggregate demands 
obtained via oracle access. 

An issue that plays a central role in our algorithms for spending constraint utilities is precision of 
prices and demands. This seems to have been treated only in minor detail in some of the previous works. 
For the algorithm of |24j in linear exchange markets these issues are discussed in depth. However, their 
solution is to change the agents side and alter utility values for maintaining bounded precision throughout 
the algorithm. However, computing an approximate equilibrium in unknown markets with demand oracle 
access seems impossible via this route. 

As a consequence, unlike for the existing algorithms in the linear case, the surpluses encountered 
by our algorithms might now become negative. Hence, additional events have to be taken into account 
upon increasing prices. Moreover, a significant challenge lies in maintaining the precision of prices to 
be polynomial throughout the algorithm. To overcome this problem, we make use of a novel tool that 
we term ratio graph. This graph is defined for a vector of prices p. The goods are the vertices, and 
we draw an undirected edge between goods j and k if the ratio of prices Pj/pk can be expressed by two 
L-bit numbers. For an intuition, observe that if some agent i has the same bang-per-buck for two goods 
j and k, then Uij/pj = Uik/pk or Pj/pk = Uij/uik, i.e., the ratio of prices can be expressed by two L-bit 
numbers. Maybe surprisingly, the broad structure of the ratio graph indeed contains enough information 
to implement algorithms for finding approximate equilibria in unknown spending constraint markets. 

1.3 Related Work 

The problem of computing market equilibria has been intensively studied, and the literature is too vast 
to survey here. We provide an overview of the work most directly relevant to ours. There is a large body 
of work on algorithms for computing equilibrium using full market information. The first combinatorial 
polynomial-time algorithm for linear Fisher markets was given by m- Later, [32| provided a polynomial¬ 
time algorithm for Fisher markets with spending constraint utilities by extending combinatorial techniques 
of |2I]. Strongly polynomial-time algorithms are also known, for Fisher markets with linear |361I44| and 
spending constraint utilities |44| . 

For linear exchange markets, [32] and [l6] obtained polynomial-time algorithms based on ellipsoid and 
interior point methods on a convex program, respectively. [ 23 | gave the first combinatorial polynomial-time 
algorithm for this problem, which was recently improved by [23|. For exchange markets with spending 
constraint utilities, [22] gave an algorithm to compute a (1 -|- e)-approximate equilibrium, for which the 
running time dependence on e is 0(l/e^). [26] gave a strongly polynomial-time algorithm for markets 
with Cobb-Douglas utilities. 

For the general case of WGS markets with unique and continuous demands, a polynomial-time al¬ 
gorithm was obtained by M- Note that this algorithm relies heavily on the ellipsoid method. For the 
Fisher setting, the famous Eisenberg-Gale convex program |28j captures market equilibrium under linear 
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utilities. m generalized it to work for any homogeneous utility functions, many of which satisfy the WGS 
property. 

As an alternative approach to compute market equilibrium, tdtonnement was defined by [IS] - al¬ 
gorithms which have access to endowment and utilities only via demand oracles. Usually, tatonnement 
procedures conduct price updates separately for each good, sometimes based on derivatives of the demand 
as a function of price. In the computer science literature, m gave a discrete tatonnement process that 
converges to an (1 + e)-approximate equilibrium for WGS markets. It has a convergence time polynomial 
in the input size and 1/e, and it does not query the original market since one needs to add auxiliary 
agents. By taking a different approach, our algorithm for the same market setting improves this rate 
to polynomial in the input size and log(l/e). Also, we only rely on approximate demand queries to the 
original market. 

More recently, m established the first fast converging discrete version of tatonnement for WGS 
markets. The convergence time depends on various market parameters. It requires a non-zero amount 
of money in the market, so it works for the special case of Fisher markets and beyond, but it is not 
applicable to the full range of exchange markets. For Fisher markets, many additional results [3lBl[9l-[TT] 
on the convergence of tatonnement processes beyond WGS markets were derived. 

2 WGS Exchange Markets 

In this section we describe the algorithm for WGS exchange markets. As previously mentioned, we 
assume that we are only granted access to an approximate demand oracle and are restricted to finite 
precision arithmetic computations. In order to make our algorithm and its analysis more accessible, we 
will simplify the problem in the remainder of this section by assuming that we are equipped with (1) 
exact real arithmetic, and (2) an exact demand oracle. This significantly simplifies the analysis in terms 
of notation and calculations, and as such concentrates on the key ideas of the algorithm. A complete and 
rigorous version of this section, presenting the entire algorithm and its proof with approximate precision, 
can be found in Appendix lAl 

We apply algorithm Alg-WGS-Precise given below. The main idea is to repeatedly identify a 
subset of goods Gi by finding a gap in the sorted order of surpluses. It then raises the prices of Gi by 
a common factor x until the surplus gap is almost closed, or the smallest surplus of Gi becomes (close 
to) 0. More formally, given price vector p = (pi,... ,Pm), value x € M+ and subset S' C G, the algorithm 
uses Update(p, x, S), which is the price vector p' = (pi,... ,p'm) with p[ = x ■ pi i G S and p( = p* 
otherwise. To implement this process, the algorithm relies on two parameters Di and D 2 based on the 
following assumptions. 

Assumption 2.1 (Bounded Price). There exists a market equilibrium (p*,x*) with 1 < p* < 2^^, Vi E G. 

Assumption 2.2 (Continuity). For any priee veetor p sueh that 1 < Pi < 2^^ for eaeh i, < 2 ^^ 
for every i,j, where Si is the surplus money of good i in 0{p), and D 2 is a polynomial of the input size. 

These two assumptions are precisely the ones from [13 ED. Assumption 12.11 about bounded prices 
is fairly mild and in many cases necessary for an efficient algorithm to compute a (strong) approximate 
market equilibrium. Assumption 12.21 about continuity is also satisfied by many natural markets, for 
example, markets with CES utilities with 0 < p < 1. Note, however, that it is not satisfied for linear and 
spending constraint markets, and hence we must develop new tools and procednres in Section |3 

To measure progress towards equilibrium we use a potential funetion <h(pt) = ||0(pt)|p. We start by 
proving a number of claims about the price vector pt. The first claim shows that with respect to exact 
demands, our algorithm monotonically reduces the 1-norm of the surpluses of all goods from 2m. 
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Algorithm 1. Alg-WGS-Precise 
Input : number of goods m, demand oracle O, precision bound e > 0 
Output: Prices p of an (1 + e)-approximate market equilibrium 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 


Set initial price po (1, 1, • • • , 1) and round index t •(— 0. 
repeat //round t 
t i — t \ 

S — ('Si) ■ ■ ■ ) 'Sm) ^ C^(Pt—l) 

Sort s such that Si^ > Si^ > ■ ■ ■ > Si^. 

Find smallest k, such that < 0 or Si^. > (1 + 

Set Gi ■(— {ii,..., ifc} and G 2 <— G\Gi 

Find the largest x, such that in s' = 0(UPDATE(pt_i, x, Gi)) it holds 
min{s' I i E Gi} = max{{s( | i E G 2 } U {0}}. 
pt ^ UPDATE(pt_i,X,Gi) 
until ||G(pt)||2 < £2 

return pt 


Claim 2.1. In Alg-WGS-Precise, |G(po)| < 2m and |0(pt)| is non-increasing in t. 

Proof. Let di be the exact demand for good i under price po, then |G(po)| = di+m = 2m. 

Next, by the criteria to define Gi and G 2 in each round, we have {i \ 0{pt-i)i < 0} C G 2 . 

During round t, only prices of goods in Gi are increased. By the WGS property, we know 0{pt)i > 
0{pt-i)i for every i E G 2 . Further, note that min{0(pt)j | i E Gi} > 0 since min{0(pt)j | i E Gi} > 0. 
Hence, we do not introduce any new negative surplus in 0{pt). Thus, we have 

\0{pt-i)\ = -2 0{pt-i)^>-2 Y. G(Pt-i)*>-2 E 0{pt)i = \0{pt)\ . 

O(pt-i)i<0 O(pt)i<0 Ci(pt)i<0 


□ 

The next two claims bound the range of prices we encounter, which is important for showing that we 
approach the unique market equilibrium. 

Claim 2.2. Throughout the run of Alg-WGS-Pregise, every good with negative surplus has price 1. 
Hence, there will be at least one good whose price remains 1. 

Proof. Observe the following three simple facts about the surplus 0{pt) resulting from exact demands: 
(1) Throughout the algorithm we never increase the price of any good with negative surplus. (2) The 
surplus of any good does not change from non-negative to negative. (3) For any non-equilibrium price 
vector, there will always be a good with negative surplus. These facts are direct consequences of the 
conditions used to classify goods based on 0{pt) in the algorithm. Together they prove the claim. □ 

Claim 2.3. In Alg-WGS-Precise, for any t>0, all prices in p^ are bounded by 2^L 

Proof. Let p* be equilibrium prices according to Assumption 12.11 We show that for any t > 0, pt is always 
pointwise smaller than p*. Assume that this is not true, let t be the smallest value such that there exists 
(Pt)i > Pi for some i. Note that according to the algorithm, we have pt = UPDATE(pt_i, x, Gi) for some 
X > 1 and Gi C [m]. Further, by our classification based on O, it is easy to see that 0(pt)i > 0 for any 
z E Gi. This means from pt-i to pi, only prices of goods in Gi are increased. Let S' = {i | (pt) * > Pi}: 
then we have S' C Gi. 
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Next, we apply a sequence of price changes to p^. First, for every j ^ S we increase (pt)j to p*. Let 
p be the new price vector and consider the surplus 0(p) resulting from exact demands. By the WGS 
property of the market, the surplus of any good in S will not decrease, hence we still have 0(p)j > 0 for 
every i £ S. The sum of all surpluses in an exchange market is always 0, so X) j^s<^ip)j < 0 . 

Now we decrease the price of every good i £ S from pj to p*. Then p becomes exactly p*. This 
process will not increase surplus of any good j ^ S. Thus we still have 0. This contradicts 

the assumption that p* are prices of a market equilibrium. □ 

The following claim establishes a relation between the surplus of a good with respect to exact demands 
before and after a multiplicative price update step. 

Claim 2.4. For any price vector p, x > 1 and S C [m], 0(Update(p, x, 5))i < x • 0(p)i for any i £ S. 

Proof. Because prices are scalable in exchange markets, we have 0{x • p) = x • 0(p) for any value x > 0. 
Also, by the WGS property, when we decrease any set of prices, this will not increase the demand for goods 
with untouched price. Since these goods have untouched price and non-decreasing demand, they also enjoy 
non-decreasing surplus. Therefore for any i £ S,we have C)(Update(p, x, S))i < 0{x-p)i = x-C>(p)j. □ 

The next lemma is the key step in the proof of our main result. It establishes a multiplicative decrease 
of the potential function at the end of many of the rounds. Let i? be a sufficiently large constant to be 
explained in the end of the proof of the following lemma. 

Lemma 2.1. Ifx< 1+;^^ the end of round t in Alg-WGS-Precise, then 4’(pt) < <k(pt_i) 

Proof. We let s = and s' = 0{pt) throughout the proof. An intuition of the proof is as follows. 

In Alg-WGS-Precise, by the conditions used to define Gi and G 2 , we always have Si^, > Si^^je and 
^ik~^ik+i > 'S 4 /("^ + l) > Sn/e(PT' + l)- Hence, roughly speaking, every good in Gi has reasonably large 
surplus, and there is a reasonably large gap between the surpluses in Gi and G 2 . Next, at the end of the 
current round, we decreased the minimum surplus of a good in Gi to either min{s^ | i £ Gi} = 0 (Case 
(1) below) or min{s' | i £ Gi} = max{s( | i £ G 2 } (Case (2) below). In both cases, the total value of 
must decrease by a factor of 1 — Q.{l/m^). 

More formally, if the algorithm proceeds to round t, then ||s|| > e^. By the definition of set Gi, we have 
Sii < (1 + ^)si 2 < • • • < (1 + < e • S 4 . Hence, > (sq/e)^ > <^{pt-i)/{me‘^) > {e/{y/me)f, 

so the surpluses of goods in Gi are similar up to a factor of e and bounded from below. Also we have 
(sife - Sife+i)^ > {sije{m + 1))2 > 4>(pt_i)/(e2(m + ifm). 

Next we relate the surplus in the beginning and the end of a round as follows. For every i £ Gi, 
by Claim 12.41 the surplus from exact demands satisfies s' < x ■ Si. Since x < 1 -|- , it holds that 

s( < (1 + We do not touch the price of any good j £ G 2 , so the WGS property implies for exact 

demands s'- > Sj. 

Now, in order to bound the change of ‘h(pt), we consider s' according to Gi and G 2 . We distinguish 
two cases. 

Case 1: max{s( | i £ G 2 } < 0 In this case the algorithm has decreased the surplus of some good in 
Gi to 0. This decrease alone brings down the potential function 4* by a factor of 1 — H(l/m). All other 
surpluses will cause an increase by a factor of at most 1 -|- 0 (l/m^). 

More formally, the contribution of goods of Gi to 4>(pt) can be upper bounded by 






Furthermore, for every i G G 2 , by the WGS property of the market, we know Sj < < 0. Thus the 

contribution of goods of G 2 to ^(pt) can be upper bounded by 


E < E 

j=k+l 


si < 


Combining the two parts 


^(Pi) = 

i=i 


E (1 


Rm^ 


s'l < 


5 ^ + 




Rm^ / 


= 1 + 


Rm^ 


($(pt-i) - s-,) 


< + 
< (1- 


1 


Rm^ 

1 

2e^m 


1 “ ^ ) ^(P‘-i) 
e^m / 


^(pt-i) 


where the last inequality holds for any m > 2 with sufficiently large constant R. 


Case 2: max{s' | i G G 2 } >0 In this case the gap between surpluses in Gi and G 2 decreases to 0. 
Below we show that the closing of this gap yields a decrease of the potential function <I> by a factor of 
1 — All other surpluses will increase by a factor of at most 1 + 0(l/m^). In combination, it 

turns out that <I> will decrease by a factor of 1 — n(I/m^). 

More formally, in this case min{s'j \ i G Gi} = max{s'- | j G G 2 }. Let = min{s'j | i G Gi} and 
SG 2 = max{s'- I j G G 2 }. For every i G Gi, let s'j = x'si — 5i where x' = (1 + and for every j G G 2 , 

let s'j = Sj + Sj. Hence 5i, 5j > 0 for all i, j. Further, we have YllLi Si = Jlih s[ = 0, hence 


and 


Now we have 


E = E - E 

2GGi i^Gi iG.G\ 

> E + E = E + E 

isGi jeG2 ieGi ieG2 

= E'^.- 

j&G2 

iSGi 


Hpt) = E = E ~ + E 




teGa 


= Y 2^'^'5* + E + E + ^j) - Y - Y + E ^3 


GeGi 


teGa 


\j&G2 


ieGi 


ieGi 


teGa 


< x'^$(pt_i) + SGa E E E + '^*'^+1 E ^3 


teGa 


ieGi 


ieGi 


J6G2 


< X 


^$(pt_l) - (si, - Y 


i^G\ 


( 1 ) 

( 2 ) 

(3) 
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< 


1 


(4) 


a;'2$(pt_i) - -(s,, - 

, 2 1 
< ( 1 + T^ + 


R^rrfi 2e^(m + l)^TO 


4’(Pt-i) 


= (1-^1;;^)) ^(P*-i 


(5) 

( 6 ) 


Here m can be derived by expanding the quadratic formula and appropriately reorganizing the terms. 
For the step from ([T]) to ([2D, in the first bracket we overestimate the quadratic terms of s into 
In the second bracket, the terms are bounded correctly using sqi for all i £ Gi and sgj for all j € G 2 . For 
the final two terms in ([ID and ([2D we use the definition of and and the fact that x' > 1. For the 
step from ([2D to ([3D , for the second bracket of ([2D we note sgi > sg 2 and the two sets Gi and G 2 have the 
same sums of 5-terms. By the same argument, we can transform the last two terms of ([2D as shown. From 
([HD to ([ID and then to ([5D, we use the bound for 5i and (s^ — > <h(pj_i)/(e^(m -|- l)^m) as 

shown above. 

Finally, using a sufficiently large constant R, the multiplicative term in ([5D can be decreased to 
strictly less than 1 for every m > 2. The hnal expression in ([6D captures the asymptotics and proves the 
lemma. □ 


Observe that the previous lemma shows a decrease in the potential only for rounds in which the factor 
X is rather small. The next lemma shows that there can be only a limited number of rounds with a larger 
value of X. 

Lemma 2.2. During a run of Alg-WGS-Precise, there ean be only 0(m'^Di) many rounds that end 
withx>l + -^. 

Proof. By Claim [231 every price can be increased by a factor of 1 -|- at most 0(log;^_,_]^/^^3 2^^) = 
0{m^Di) times. Hence there can be at most 0{m^Di) many rounds with x > 1 -|- Q 

Finally, we can assemble the properties to show the number of rounds to reach an (1 -|-e)-approximate 
market equilibrium is polynomially bounded. 

Lemma 2.3. For any market that satisfies Assumptions \2A\ and \2.2l and for any e > 0, Alg-WGS- 
Precise returns the priee veetor of an (1 -|- e)-approximate market equilibrium in a number of rounds 
polynomial in the input size and log(l/e). 

Proof. Let xt be the value of x we find in round t of Alg-WGS-Pregise. First because at least one 
price will increase by a factor of xt in round t, by Claim 12.31 we have Ht < 2™^!. At the end of 
round t, if x* > 1 -|- let s = max{sj | Si € 0{pt-i)} and s' = maxjsi | Si S 0{pt)}, then by 

Claim [2~i] we have <h(pt) < ms'‘^ < mx‘fs‘^ < mxj‘h(pt_i). Moreover, by Lemma [2.21 there will be at most 
0{m‘^Di) such rounds. Hence the total increase of <h(pt) in these rounds will be no more than a factor 
of nx,>l+l/Rm3 mx? < ^ ^Oim^D,), 

For all other rounds, we have x < 1 -|- and by Lemma (231 the potential function is decreased by 
a factor of 1/(1 — H(^)). Therefore the total number of rounds before <h(pt) < will be at most 

/ Tjf>{m^Dx)\ / 3 1 

O ( - 1 = O iDimAogm + m^log- 

□ 
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The final step is to argue that the total running time (in terms of oracle queries) of the algorithm is 
polynomial. Given Lemma [2.31 it remains to show that each round only invokes polynomially many oracle 
queries. Intuitively, this can be achieved by finding x in line 8 of Alg-WGS-Precise via binary search. 
However, for a formal proof it is unavoidable to make statements about the required bit precision of p, 
X, and the approximate demand oracle, which we did not discuss here. In the following we present the 
statement of the main theorem. The formal proof, together with the statement of algorithm Alg-WGS 
based on bounded precision and full details on its analysis, is deferred to Appendix [Al 

Theorem 2.1. For any market that satisfies Assumptions \2A\ and \2.Sl and for any e > 0, Alg-WGS 
returns the priee veetor of an (l + e)-approximate market equilibrium in time polynomial in the input size 
and log(l/e). 

Remark 2.1. Our main goal in the analysis was to establish a bound on the running time that is 
polynomial in m, L, and log(l/e). For the sake of simplicity we did not optimize the bounds beyond 
being polynomial. It appears that the dependence on m can be significantly improved, e.g., by a more 
precise analysis of the actual number of rounds with a: > 1 -|- and their impact on the potential. 
Moreover, based on our preliminary experiments, it appears that using pL = 0(e/m^) is sufficient, and 
the algorithm converges to an equilibrium much faster than the bound predicts. 

3 Exchange Markets with Spending Constraint Utilities 

In this section we discuss our algorithm for exchange markets with spending constraint utilities. Spending 
constraint utilities are defined in |221I42| . where the utility derived by agent i from good j is given by a 
piecewise linear concave (PLG) function fij. The overall utility of agent i is additively separable among 
goods, i.e., tti(x) = J2j£G Each fij is a PLG function with a number of linear segments. Each 

segment k has two parameters: the rate of utility Uijk per unit of good derived on segment k and the 
maximum fraction Bijk of budget that can be spent on segment k. All Bijj. are strictly positive, and 
concavity implies Uijk > Here we assume all UijkS are integers, all Bijk,Wij^s are rational 

numbers and the whole input can be represented in no more than L bits. 

Spending constraint markets may not have an equilibrium m, however under mild conditions, there 
is always a rational equilibrium [35]. Henceforth we will assume the following sufficient condition. Let 

r(5) = {j e G \ Wij > o,i e S}. 

Assumption 3.1 (Sufficiency Gondition). For any subset S of agents, ifV{S) G then there exists 
i & S and j & G \ r(5) sueh that Uiji > 0. 

Let us first characterize the demand of each agent i under spending constraint utilities. Given nonzero 
prices p, define the bang-per-buck relative to p for segment k in to be Uijk/pj- Sort all segments of 
agent i by decreasing bang-per-buck value, and partition them by equality into classes: Qi, Q 2 , ■ ■ ■, Qt- 
Then an allocation is a demand bundle of agent i if and only if there is an integer ti such that all 
segments in partitions Qi,..., Qp-i are all fully allocated, and no segments in partitions Qti+i,Qti+ 2 , ■ ■ ■ 
are allocated. Furthermore, the total money spent on partitions Qi,... ,Qp-i is no more than agent 
z’s total budget ruj = term Qt- agent z’s eurrent partition, and Qi ,..., Qp-i agent Vs 

allocated partition. Let spentf denote the total money of agent i spent on allocated segments, and let 
spent^ denote the total money spent on allocated segments of good j. Agent i can freely demand any 
segments in her current partition, fully or partially, until her remaining budget mt — spentf is exhausted. 

Equality and Ratio Graph The main tool for the analysis of previous algorithms in spending con¬ 
straint markets is an equality graph, denoted by EG{'p). This graph remains completely unknown to 
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our algorithm, but it proves useful when proving properties of the convergence process. The vertex set 
of this bipartite graph consists of the set of agents A and the set of goods G. Given a price vector pt, 
we introduce an edge from agent i to good j if and only if agent i’s current partition Qf- contains one 
segment that belongs to utility function fij for good j. The edges in EG(p) are called equality edges. 
Observe that this graph changes throughout the process when we update the price vector p*. 

Based on ^^^(p) we can construct an equality network denoted by iV(p): First for each edge in 
EG{p) from agent i to good j, let k be the corresponding segment for fij that belongs to Qt^. We assign 
a capacity of Cij = Bij^rui to this edge. Next add a source vertex s and a sink vertex t. For each agent i, 
add an edge from s to z with capacity m* — spenti, and finally add an edge from every good to sink t with 
infinite capacity. It is easy to see that every maximum flow in A^(p) corresponds to a feasible demand 
allocation for each agent. 

Similar to [2T], we define a balanced flow as a maximum flow in A^(p) that minimizes + 

spent^j — Pj)'^, where ijt is the flow along edge (j, t), which also denotes the money spent on good j on 
segments of current partitions, and spent^j is the amount of money spent on allocated partitions of good j. 
By assumption, 0{p) returns the surplus vector derived from (any) balanced flow of the network A^(p)ll. 

As mentioned above, our algorithm cannot see EG{p). We use a different structure that can be 
observed on the goods side. 

Definition 3.1. The ratio graph RG{M,p) is an undirected graph with m vertices (where m is the 
number of goods in the market), and for any two goods i and j, {i,j) is an edge if and only if Pi/pj can 
be represented as a ratio of two integers, each of value at most 2^ — 1. 

We can compute the ratio graph using only the price vector and the input size bound. It allows us to 
retrieve some information about the hidden structure of EG{p). 

Claim 3.1. Let L be the upper bound on the number of bits to represent each utility parameter and p 
be a price vector. For any price vector p and goods i,j that are connected in EG{p), i and j are also 
connected in RG{M, p) for any M > L. 

Proof. If good i and good j are connected in EG{p), then there exist goods i = io,ii,... ,ik-iAk = T 
such that for each t < k, there exists some agent at that has the same bang-per-buck for two segments 
for goods it and it+i, respectively. Then we have PitIPit+i = Uatitki/uatit+ik 2 some ki and /c 2 . This is 
a ratio of two integers, each of value at most 2^ — 1. This implies {it,it+i) G RG{M,p), so i and j are 
connected in 1?G(M, p). □ 

There exists an algorithm for computing exact equilibrium prices in Fisher markets with spending 
constraint utilities [52]. [22| give an FPTAS for exchange markets with spending constraint utilities. The 
algorithm finds an (1 + e)-approximate equilibrium in time polynomial in the input size and 1/e. 

For spending constraint markets, we extend our approach for WGS markets. The challenge is that 
surplus can change in a non-continuous way when prices change the current partitions of the agents. 
However, we show how to use the linear structure of the market to get rough information about these 
breakpoints. Also, we maintain prices within a polynomial precision and guarantee convergence to an 
approximate market equilibrium. Finally, when 5>(p) becomes small enough, we convert the approximate 
equilibrium to an exact one using a procedure Alg-Spending-Exact. 

Our only assumption is that the whole input can be represented within L bits, and L is known 
to the algorithm. This implies as a corollary a variant of Assumption 12.11 “ there is an exact market 
equilibrium with prices p s.t. < 2^^ where Di is a polynomial of m and L. As mentioned before, 

mini Pi 

Assumption 12.21 does not hold in spending constraint markets. 

®It can be easily shown that every balanced flow gives the same surplus vector. 
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Algorithm 2. Alg-Spending: Framework for Spending Constraint Markets 
Input : number of goods m, demand oracle O, number of bits L to represent the whole input 
(including each Uijk,Bijk,Wij)\ solution precision e 
Output: Prices p of a (1 + e)-approximate market equilibrium 
Parameters: e' = 77 - 1 = 

2y/m 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 


Set initial price po ( 1 ,1, • • • , 1) and round index t •(— 0 . 
repeat //round t 
t i — t \ 

S — ('SI5 ■ ■ ■ ) Sfn) ^ C^(Pt—1) 

Sort s such that sq > Sjj > • • • > • 

Find smallest k, such that 84,^^ < 0 or Si^. > (1 + 

Set Gi •(— {ii,... ,ifc} and G2 [m]\Gi 

Binary search the smallest x G (l,oo), such that in s' = 0 (UPDATE(pt_i, x, Gi)) it holds 
min{s' I i € Gi} < max{{s' | i € G2} U { 0 }}. 

Pt ^ UPDATE(pt_i,X,Gi) 
until ||G(pt)||2 < e'2 

return p^ 


3.1 The Framework 

Alg-Spending specihes the general framework of our algorithm, which is similar to the approach taken 
previously for WGS markets. As mentioned in Section 11.11 here we do not resort to approximation 
parameter /i, but instead compute the exact surplus. We first analyze Alg-Spending and show that 
it needs only a polynomial number of rounds to converge to an approximate market equilibrium. For 
now, our analysis disregards all precision and representation issues. In particular, we assume to find the 
exact value x using binary search, irrespective of the number of bits needed for representation. Also, the 
update of prices from pt_i to pt will multiply all prices of goods in Gi by x, irrespective of the number 
of bits required to represent them. In our final algorithm below, we will show how to address these issues 
to obtain a (true) polynomial-time algorithm. 

The analysis of Alg-Spending proceeds roughly as in the previous section. We rely on the following 
lemma. 

Lemma 3.1. The demands returned by G(p) satisfy the WGS property. 

Proof. We make use of a max-min fair property for balanced flows in linear markets proved in previous 
work. A vector s is called max-min fair iff for every feasible vector s' and i such that Si < s[, there is 
some j with Sj < Si, such that Sj > s'-. The following claim is proved by [21| . 

Claim 3.2. \21^ The surplus veetor of a balaneed flow in N (G) is max-min fair among all feasible surplus 
veetors. 

Although this claim is for linear markets, it can also be directly applied to spending constraint markets 
because the network flows are designed only with respect to the current partition of each agent i. Within 
this domain the spending constraint market behaves exactly like a linear market. 

We proceed to prove Lemma 13.11 by contradiction to this claim. Suppose we increase the price of 
some good k from pk to pk -\- 5. We denote the old prices by p and the new prices by p'. Let s = 0{p) 
and s' = O(p'). Now assume for contradiction that there exists some i f k such that < si. Let 
5 'g = {j £ G* I Sj > S£,j f fc} and 5^ = {* G A | there exists j G Sq such that fij > 0 in A(p)}. It is 
easy to verify the following properties. 
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(a) For any edge {i,j) in A^(p) with j ^ k, let seg be the segment in agent i’s current partition that 
corresponds to this edge. Then with new prices p', either seg belongs to the allocated partition 
(i.e., is fully allocated), or seg is still in agent i’s current partition (i.e., (i,j) also exists in A^(p')). 

(b) For any edge (i,j) € iV(p') with i G Sa and j 0 Sq, let seg be the segment in agent i’s current 
partition that corresponds to this edge with prices ph If seg is also in i’s current partition with 
prices p, then seg is fully allocated with price p (otherwise agent i can reroute some flow from Sq 
to this segment to obtain a more balanced flow in A^(p)). If seg is not in i’s current partition with 
price p, then all segments in i’s current partition with price p must be fully allocated in the demand 
allocation with prices ph 

The above two observations imply that with prices p', we can rearrange flows from Sa to Sq to get a 
new feasible surplus vector s" such that for all j & Sq, Sj > Sj, and for all other edges j ^ Sq, s” = s'-. 
In particular, we have s'- > Sj > s^ > s^ for all j E Sq- This contradicts the fact that s' is max-min fair 
and proves Lemma l3. II □ 

Now the following properties can be proved using literally the same proofs as for WGS markets before. 
Claim 3.3. For Alg-Spending the following properties hold. 

1. In Alg-Spending, |C>(po)| < 2m, and |0(pt)| is non-increasing in t. 

2. For any price vector p, x > 1, and S C [m] it holds C>(Update(p, x, S))i < x ■ 0(p)i for any i £ S. 

3. The number of rounds that end with x > 1 -|- in Alg-Spending is 0{m^Di), for a sufficiently 
large constant R. 

We also show a version of Claim E2] for spending constraint markets, which needs some extra work. 
Unlike for WGS markets, the surplus of a good can change from non-negative to negative. Thus, the 
proof of Claim [22] does not directly transfer to spending constraint markets. Instead, we first show the 
following. 

Claim 3.4. Let St be the set of goods with negative surplus and price strictly greater than 1 at the end 
of round t in Alg-Spending. For any T C St, let T{T,pt) be the neighbors of set T in EG{pt), i.e., 
r(T,pt) is the set of agents who are interested in at least one good in T under price pt. Let B{T{T,pt)) 
be the sum of budgets of agents in F(r, pt). Then we have B{T{T,pt)) > Z)ieT(Pi)i- 

Proof. We prove this claim by induction. The claim is trivially true for round 0. Assume that it is true 
for any round t <t', then at the end of round t = t' 1, consider two cases: 

• min{s' I i E Gi} > 0. Because the surplus of any good in G 2 is non-decreasing in round t, we 
have St C St-i. Further, the algorithm does not increase the price of any good in T in round t. 
Hence, we also have F(r,pt_i) C T{T,pt). By the induction assumption, we have B(T{T,pt)) > 
H(F(r,pt_i)) > EiMPt-i)i = 

• min{s( I i E Gi} < 0. Assume that during this round we start from x = 1 and increase x continu¬ 
ously until it reaches its final value. We also assume that the equality graph £'G(UPDATE(pt_i, x, Gi)) 
and the corresponding balanced flow are implicitly being maintained throughout the process. For 
each agent i and any moment during this round, let F(i) denote the neighbors of agent i in the 
equality graph. There are two cases: 
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Case (a): r(z) n Gi 7^ 0 and r(z) H G2 7^ 0 - This means agent i’s current partition contains 
segments of goods in both Gi and G 2 . When we continue to increase x (and consequently the prices 
of goods in Gi), the segments of goods in Gi will have worse bang-per-buck value than segments 
of goods in G 2 . Hence, they will be removed from agent i’s current partition. Furthermore, for any 
good ji G r(i) n Gi and j2 G r(i) n G2, because we have Sjj > Sj^ and the balanced flow condition, 
it cannot happen that > 0 and < Cij^. Otherwise agent i would be able to route some flow 
from edge (i,ji) to edge (i,j2) to reach a more balanced flow. Thus, two possibilities remain: 

— = Cij 2 for every j 2 G r(i) n G 2 . This means all segments of goods in G 2 are already fully 
allocated. When x increases, the new current partition of agent i will only contain segments 
of goods in r(i) n Gi, and the surpluses will change continuously with the change of x. 

— iij^ = 0 for every ji G r(i) n Gi. This means when x increases, the segments that are being 
removed from i’s current partition are all unallocated in the current allocation. Hence, the 
surpluses will again change continuously at the current point. 

In summary, in this case the change of the surpluses will always be continuous in the change of 
X. Note that when the surpluses are changing continuously. Case 2 cannot occur before Case 1 
happens. Therefore, the Claim follows via Case 1. 


Case (b): r(i) C Gi or r(i) C G 2 . This means we will always change the prices of all segments 
in i’s current partition by the same rate. Hence, if set Qt^ changes, it can only be merged with some 
segments in Qti-i or In either case, the final value of x must be at a point where at least one 

new edge emerges in the equality graph EG{'pt) from r(Gi) C A (the set of agents incident to at 
least one good in Gi) to G 2 . Suppose we alter the equality graph by removing these emerging edges 
from HG(pi). If we recompute the balanced flow, then min{s( | i G Gi} > max{{s( | i G G2}U{0}}. 
For any T C let Ti = Tn Gi and r 2 = Tn G 2 . In this new graph, let F'(ri) be the set of agents 
who have positive flow to at least one good in Ti, and let F'(r 2 ) be the set of agents incident to at 
least one good in T 2 . Since min{s' | i G Ti} > max{s( | i G r 2 }, by the balanced flow condition 
we know F'(ri) n F'(r 2 ) = 0. Also, we have H(F'(ri)) > J2ieTiiPt)i because every good in Ti has 
positive surplus, and H(F'(r 2 )) > X]ier 2 (Pt)* because the claim is true in round t — 1. Combining 
these two inequalities gives us i?(F(r, pt)) > 

□ 

Corollary 3.1. At the end of eaeh round t in Alg-Spending, for any good i with negative surplus and 
priee greater than 1, there exists another good j with priee 1 that is eonneeted to i in EG{'pt) 

Proof. Assume for contradiction that the statement is false. Let T be the set of goods with negative 
surplus and connected with good i in EG{pt). By the balanced flow condition, none of the agents in 
F(r, pt) can have any positive flow to goods outside set T. Thus we have 0 > — -®(b(T', Pt)) — 

(pt)j. This contradicts Claim ElU □ 

We obtain the following corollary, an analog of Claim 12.21 for spending constraint markets. 

Corollary 3.2. Throughout the run of Alg-Spending, there will be at least one good whose priee 
remains 1. 

The set of properties shown so far allows to establish the following lemma, which is the key step 
for observing convergence to equilibrium. It can be seen as an adjustment of Lemma 12.11 to spending 
constraint markets. As before, it involves a sufficiently large constant R. 
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Algorithm 3. RouNDiNG (p,M): Rounding procedure 
Input : price vector p in which minj pj = 1, rounding bound M > L 
Output: Rounded price vector p' 


1 

2 

3 

4 

5 

6 

7 

8 


Let R = {| I o, 6 S ,a, b < 2^} and p' ■(— p. 
while RG{M,p') is not connected do 

Let Cl, C 2 ,..., Cfc be the connected components of RG{M, p') 

Assume without loss of generality that Ci is a component with m\ni^C\r\AP'i = 1 
For every i,j, let rjj = p'-jp'j. 

Let bij = min{^ \ x ^V,x >rij} and Bij = min{6j'j' | i' € Ci,j' € Gj}. 

Let {i,j) E argminj>j Bij. 
p' ^ Update(p', Bij,Gi) 


9 return p' 


Lemma 3.2. Ifx< (1+-^^) o,t the end of round t in Alg-Spending, then ‘h(pt) < <h(pi_i) — R (;^)) • 
Proof. Our proof uses the arguments of the proof for Lemma l2.ll We consider two cases: 

Case 1: min{s' | i E Gi} = max{{s'j | i E G 2 } U {0}} This case can be verified by observing that 
Claim IAAI holds with p, = 0. Then the proof follows using exactly the same proof as for Lemma l2.ll 

Case 2: min{s' | i E Gi} < max{{s' | i E G 2 } U {0}} Using the same argument as in the proof of 
Claim 13.41 one can show that in this case x must be at a point where at least one new edge emerges 
in the equality graph EG{pt) from r(Gi) C A. Without these emerging edges, we have min{s^ | i E 
Gi} > max{s^ | i E G 2 }. This means that we can further reduce the flows along the new edges to get 
another feasible flow in N(pt), such that with the resulting surplus vector s" we have min{sf | i E Gi} = 
max{s" I i E G 2 }. 

Next, using again the same proof as for Lemma l2. 11 for a sufficiently large constant R it follows that 
||s"||2 < ^>(pi_i) (1 - 0 (;^)). Hence, $(pt) < ||s"f < ^(pt_i) (l -H □ 

3.2 Precision and Representation 

It is now tempting to think that using a similar argument as in the general WGS case, Claim 13.31 
Corollary 13.21 and Lemma [3.21 provide an algorithm that converges to an approximate market equilibrium 
in polynomial time. However, an issue arises with regard to the precision and representation of the 
prices: In each round, x could potentially be a rational number involving prices and surpluses, and after 
multiplying each price in Gi by x, the bit length to represent a price can double in one round. This means 
that after a polynomial number of rounds, we may require an exponential number of bits to represent the 
prices of some goods as well as the desired factor x. 

Recall the ratio graph and Claim [3Tl As an adjustment, we run the following procedure Rounding(p, M) 
at the end of each iteration of the main loop in Alg-Spending. Its purpose is to round the prices within 
polynomial bit length while maintaining the structure of equality and ratio graphs. Thereby, we will show 
that the value of potential function 4>(p) will not be increased dramatically. 

Our new algorithm Alg-Spending-Rounding is simply the framework Alg-Spending with the 
following modifications: 
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(1) Set M = log 2 In each round, binary search x within domain V = {f \ a > b,a,b ^ , a, b < 

22 mM+Lj^ instead of (l,oo). 

(2) At the end of each round t, update pt ^ RouNDiNG(pt, M) with M as above. 

Lemma 3.3. Given any price vector p, Rounding(p, M) terminates in time polynomial in m, M. The 
returned price vector p' satisfies: 

(a) there exists i ^ G with pi = p[ = 1, 

(b) every price can be represented as a ratio of two integers, each of value at most 2 ™^, 

(c) EG{p') contains every edge present in EG{p), and 

(d) Pi < p'i < Pi + for every i ^ G. 

Proof. Property (a) holds since we never increase the price of goods in set Ci, and property (b) can be 
derived based on property (a) and the fact that RG{M, p') is connected when the algorithm terminates. 

Next we claim that for any i,j, there does not exist any x € R, such that Pi/pj < x < p'ijp'j. This 
is because by design, pi/pj has to reach x in some iteration before it grows beyond x. But starting from 
that moment until the end of the algorithm, i and j will be connected in RG{M,p). The ratios between 
two prices in the same connected component in RG{M,p) remain unchanged. Hence Pi/pj will never 
grow beyond x. 

This claim also proves property (c). By Claim IHTl a pair of goods i and j connected in EG{p) remain 
connected in RG{M,p). Hence, their ratio of prices will remain the same in p', and they are connected 
in EG{p'). 

For property (d), according to the algorithm no price will decrease from p to p^ Next, number the 
goods such that pi = p^ = 1. For any i 7 ^ 1, let Xi = min{x € R | x > pi\. Then we have p'^/pi = p'i < Xj, 
as well as Xi < pi + 2~^. Therefore pi < p\ < Pi + 2~^. 

Finally, in each iteration we add at least one edge between two connected components in RG{M,p). 
Thus the algorithm will terminate after at most m — 1 iterations, and it is easy to check that each iteration 
runs in polynomial time. This proves the claim. □ 

Once the prices are bounded by a fixed polynomial bit length, we can also bound the length needed to 
encode the desired x in each round. This implies that we can find x in the framework in polynomial time 
using binary search. 

Lemma 3.4. In every round of Alg-Spending-Rounding, the desired x can be represented as a ratio 
of two integers, each of value at most 2 "^^+^+^ log m _ ^ _ 

Proof. Let p = pt_i and consider the structure of EG{p). We let Ai = F(Gi) be the agents connected to 
goods in Gi. Also, let A 2 = A \ Ai. There is no {i,j) G EG{p) with i G Ai and j G G 2 , since otherwise 
0{pt-i) could increase the money spent on goods in G 2 and further decrease <I>(pi_i). For simplicity, we 
will also assume that there is no edge (i,j) G A 2 x Gi, since no agent spends money along these edges 
and they immediately disappear once we start increasing prices in Gi. 

Now we increase p on goods in j G Gi by x and get a new price vector p(x). This only generates new 
edges (i,j) G Ai x G 2 . Furthermore, we drop only edges {i,j) G Ai x Gi. To verify this let us consider 
the other possibilities. The relation between marginal utility values uijk/pj and Uiyi^ijpy for goods in 
G Gi does not change, since both pj and pfc are both multiplied by x. Hence, there are no new edges 
(i,f) G Ai X Gi. For the same reason, there are no new edges {i,j) G A 2 x G 2 . The bang-per-buck 
of goods in Gi decreases, so we also do not introduce edges (i,j) G A 2 x Gi. In fact, this also implies 
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that we do not drop any edges (i,j} G A 2 x G 2 - prices and bang-per-buck relations among goods in 
G 2 do not change at all and G 2 becomes more attractive compared to Gi. Finally, there exist no edges 
{i,j) € (^1 X G2) U {A2 X Gi) that could be removed. This shows that we only generate new edges 
between Ai and G 2 , and we only drop edges between Ai and Gi. 

For any given x, consider the residual graph of N (Update(p, x, Gi)). Let G be an arbitrary connected 
component in this graph. Let Ga be the set of agents in G and Gg be the set of goods in G. Then we 
know that all goods in Gg have the same surplus, and all flow going through Gg comes from agents in 
Ga- This implies the following equation: 

^ Pi + x pi= Pi + x ^ pi + IGcIs, 

i&CA\Gi ieCAnCi ieCcVCi ieCcnGi 

where s is the surplus of (any) good in this component. 

Now let us focus on the moment where x reaches the desired value at the end of round t according to 
the algorithm. At this moment, one of the following properties must hold: 

(1) min{s' I i G Gi} = 0. Then for the connected component that contains a good of surplus 0, we 
have s = 0 in the above equation. All initial prices pi are ratios of integers with values at most 
2 mM, gQ when we solve the equation for x, the solution is a ratio of two integers with value at most 

( 2 ) min{s' I i G Gi} = max{s'- | j G G 2 }. In this case, we have two possibilities: 

— There exist two connected components in the residual graph of A^(Update(p, x, Gi)) that have 
the same surplus. Applying the above equation to these two components, we can solve for x, 
and the solution will be a ratio of two integers with value at most 

— A new edge {i,j) G Ai x G 2 appears, then for agent i good j G G 2 becomes equally attractive 
as some k ^ Gp. Uij/pjX = Uik/pk, or, equivalently, x = {uijPk)/{uikPj)- Pij and pk can be 
represented as ratio of integers of value at most 2™^ by Lemma and Uij and Uik are both 
integers of value at most 2^ — 1. Hence, every value of x at which a new edge evolves in 
£'G(p(x)) can be presented as a ratio of integers of value at most 


□ 


We now bound the impact of replacing p^ by RouNDiNG(pt, M) in the function ‘h(pt). 

Lemma 3.5. <h(RouNDiNG(pt, M)) < <h(pt) + for any round t. 

Proof. Let p( = RouNDiNG(pi, M). By Lemma f3.3f c). we know this rounding procedure does not remove 
any edges in EG{pt). Let / be a balanced flow in N{pt). Then by Lemma IH^ di we can construct a 
feasible flow /' in A'(RouNDiNG(pt, M)), such that fij < fk < fij + 2~^ for every i,j. Let s be the 


surplus vector derived from f, then we have 

Si — 0{pt)i\ < m2 ^ for every i 

\Hl-\\oipt)g = 

- 0{ptfi) 

< 

J2i‘2m2-^\Oipt)i\+m^2-^^) 

= 

2m'^2-’^\0{pt)\ +m^2-^^ 

< 

W2"^ + 

< 

5m^2-^ . 
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Note that s is just one feasible surplus vector for price vector p^, and <h(Pt) minimizes the t'2-norm of 
surpluses among all feasible surplus vectors. Hence, <l>(RouNDiNG(pt, M)) < ||s||2 < \\0{pt)\\2 + 5m^2~^. 
This proves the lemma. □ 

Rounding can be used to obtain an algorithm that converges to an approximate market equilibrium 
in polynomial time. 

Lemma 3.6. For any spending constraint exchange market satisfying Assumption IRil an (1 + e)- 
approximate market equilibrium can be computed in time polynomial in m, L, and log(l/e). 

Proof. In the Alg-Spending framework with Alg-Spending-Rounding, we know by Lemma [3]2] that 
at the end of each round t and before calling Rounding, ‘h(pt) < <h(pt_i)(l-R(^)). Ifcl>(pi)>e'^ 
we have /mf < d>(pt)/m^. Thus by Lemma 13.51 

$(ROUNDING(pt,M)) < (l + ^)$(pi) 

This implies we can employ the same proof as for Theorem l2.1l to show that after finishing Alg-Spending- 
Rounding, we arrive at a (1 + e)-approximate market equilibrium. Because M is a polynomial in the 
input size and log(l/e), the binary search and the Rounding procedure run in polynomial time in each 
round of the framework. Hence, the running time is polynomial in the input size and log(l/e). □ 

Finally, it remains to convert the approximate equilibrium to an exact one. To achieve this, we rely 
on full information about the spending constraint utilities. While this step can be seen as an extension 
of the technique developed in [2Tj for the linear exchange markets, there are several challenges due to 
the much more involved setting of spending constraint utilities, where the allocated partitions make the 
remaining budgets of agents and the values of goods dependent on too many parameters. Using a more 
involved procedure we are able to handle the extra complexity of the problem. Our result resolves an 
open question of [M] of hnding an exact polynomial time algorithm for exchange markets with spending 
constraint utilities. A detailed discussion of this final step can be found in Appendix [B1 This yields the 
hnal theorem in this section. 

Theorem 3.1. For any spending constraint exchange market satisfying Assumvtion \3. 11 Alg-Spending- 
Ex AGT returns the price vector of a market equilibrium in time polynomial in m and L. 
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Algorithm 4. Alg-WGS 


Input : number of goods m, approximate demand oracle O, precision bound e > 0 
Output: Prices p of a (1 + e)-approximate market equilibrium 
Parameters: u = 

^ H\r 




^ _ j_^2^l+^2+logm'| 


£' = 


2y/rn 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


Set initial price po (1,1, • • • , 1) and round index t •(— 0. 

Let V = {^ \ a > b,a,b ^ ,a,b < A} 
repeat //round t 
t i — t \ 

s = ^0(pt_i,/i) 

Sort s such that Si^ > Si^ > ■ ■ ■ > Si^. 

Find smallest k, such that or Si^. > (1 + 

Set Gi ^ {ii,... ,ik} and G 2 G\Gi 

Binary search the largest x gV, such that in s' = C>(UPDATE(pt_i, x, Gi)) it holds 
min{s' I i E Gi} > max{{s' | i E G 2 } U {/i}}. 

Pi E- UPDATE(pi_i,X,Gi) 
until ||G(pi)|p < e'2 

return pi 


A WGS Exchange Markets 


In this section we describe the complete algorithm Alg-WGS for WGS exchange markets. Recall the 
assumptions from Section [21 which we restate here for completeness. 

Assumption [2TT] There exists a market equilibrium (p*,x*) with 1 < pt < E G. 

Assumption 12.21 For any priee veetor p sueh that 1 < Pi < 2^^ for eaeh i, 1I < 2^^ for every i,j, 

Opj 

where Si is the surplus money of good i in G(p), and D 2 is a polynomial of the input size. 

Throughout the analysis and proofs below, if s = 0{p) for some p, we use s to denote the surplus 
vector returned by the //-approximation demand oracle with the same price vector, i.e., s = G(p,//). We 
proceed along similar lines as in Section |21 and the proofs of the first claims closely resemble the versions 
for the exact oracle. For completeness, we provide them here for the approximate oracle. 

Claim A.l. In Alg-WGS, |0(po)| < 2m and |G(pi)| is non-increasing in t. 

Proof. Let dt be the exact demand for good i under price po, then |0(po)| = J2i \di — l\ < J2i di+m = 2m. 
Next, by the criteria to define Gi and G 2 in each round, we have {i \ < 0} C G 2 : To see 

this, observe that the surplus resulting from the approximate 0{pt-i) differs by at most an additive 
/i = e/{RimJ), so a good i with 0(pt-i )i<0 will always be classified in G 2 with respect to 0(pt-i)i. 

During round t, only prices of goods in Gi are increased. By the WGS property, we know 0{pt)i > 
0{pt-i)i for every i E G 2 . Further, note that min{0(pj)j | i E Gi} > 0 since min{0(pt)j | i E Gi} > h- 
Hence, we do not introduce any new negative surplus in 0{pt). Thus, we have 

|G(Pi_i)| = -2 ^ G(pt_i), >-2 ^ G(pt_i), >-2 ^ 0{pt)^ = \0{pt)\ . 

C>(pt-i)i<0 O(pt)i<0 C>(pt)i<0 


□ 


The next two claims bound the range of prices we encounter, which is important for showing that we 
approach the unique market equilibrium. 
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Claim A. 2. Throughout the run of Alg-WGS, every good with negative surplus has priee 1. Henee, 
there will be at least one good whose priee remains 1. 

Proof. Observe the following three simple facts about the surplus 0(pt) resulting from exact demands: 
(1) Throughout the algorithm we never increase the price of any good with negative surplus. (2) The 
surplus of any good does not change from non-negative to negative. (3) For any non-equilibrium price 
vector, there will always be a good with negative surplus. These facts are direct consequences of the 
conditions used to classify goods based on 0(pt) in the algorithm. Together they prove the claim. □ 

Claim A. 3. In Alg-WGS, for any t > 0, all priees in pt are bounded by 2^i. 

Proof. Let p* be equilibrium prices according to Assumption 12.11 We show that for any t > 0, pt is always 
pointwise smaller than p*. Assume that this is not true, let t be the smallest value such that there exists 
(Pt)i > Pi for some i. Note that according to the algorithm, we have pt = UPDATE(pt_i, x, Gi) for some 
X > 1 and Gi C [m\. Further, by our classification based on O, it is easy to see that 0(pt)i > 0 for any 
i G Gi. This means from pt_i to pt, only prices of goods in Gi are increased. Let S' = {i | (pt) * > Pi}: 
then we have S C Gi. 

Next, we apply a sequence of price changes to pt. First, for every j ^ S we increase {pt)j to p^. Let 
p be the new price vector and consider the surplus 0{p) resulting from exact demands. By the WGS 
property of the market, the surplus of any good in S will not decrease, hence we still have G(p)i > 0 for 
every i G S. The sum of all surpluses in an exchange market is always 0, so 

Now we decrease the price of every good i G S from pi to p*. Then p becomes exactly p*. This 
process will not increase surplus of any good j ^ S. Thus we still have 0. This contradicts 

the assumption that p* are prices of a market equilibrium. □ 

At this point, let us recall Claim [23] to establish the relation between the surplus of a good before 
and after a multiplicative price update step. It does not involve the approximate oracle. 

Claim [274] For any priee veetor p, x > 1 and S C [m], G(Update(p, x. S'))! < x • 0{p)i for any i G S. 

Next, we establish a statement about the surpluses at the end of each round, which was not necessary 
for the version with exact oracles and precision. Intuitively, we increase the prices of goods in Gi until 
the minimum surplus in Gi reaches the maximum surplus in G 2 or 0. Note that p, is very small and can 
be thought of as 0. The main complication here is that we need to work with /r-approximation demands 
in the algorithm and the resulting surpluses s'. 

Claim A. 4. At the end of eaeh round in Alg-WGS, min{s' | i G Gi} < max{{s' | i G G 2 } U {^}} -|- Qp. 

Proof. According to the binary search procedure, we know that if we increase prices in Gi by a factor of 
X, then s' satisfies the condition min{s' | i G Gi} > max{{s' | i G G 2 } U {//}}. Furthermore, an increase 
by x"'' = min{y GV\y>x}<x+-^ would result in a surplus vector that does not satisfy this condition. 
Let s+ = 0(UPDATE(pt_i, x+, Gi)). By Assumption 12.21 we have 

2 ^ 2 1 1 1 2^2 - 1-^1 -i-iog m 

|s+ - s'i\ < |s+ - s'l +2p< 2^^ • (x+ - x)|pt_i| +2p< - + 2p< -- + 2p = 3p 

for every i, where the last inequality is derived by Claim ESI Thus 

min{s' I i G Gi} < min{s^ | i G Gi} + 3p < max{{s+ | i G G 2 } U {/i}} -|- 3p 

< max{{s' I i G G 2 } U {//}} -|- 6 p 


□ 
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We are now ready for the key lemma in the proof of the main result - the multiplicative decrease of 
the potential function at the end of each round. Let R 2 be a sufficiently large constant to be explained 
in the end of the proof of the following lemma. 

Lemma A.l. If x < 1 + the end of round t in Alg-WGS, then ‘h(pt) < <h(pt_i) ^1 — 17 

Proof. We use the following notation. Let s = s = and s' = 0{pt), s' = 0{pt,h)- 

The intuition of the proof is similar to the version with exact precision. By the conditions used to define 
Gi and G 2 , we always have > Sj^/e and > Sif,/{m + 1) > Si^fe{m + 1). Hence, roughly 

speaking, every good in Gi has reasonably large surplus, and there is a reasonably large gap between the 
surpluses in Gi and G 2 . Next, at the end of the current round, we decreased the minimum surplus of a 
good in Gi to either min{s^ | i E Gi} « (Case (1) below) or min{s' | i E Gi} max{s( | i E G 2 } (Case 
(2) below). In both cases, the total value of must decrease by a factor of 1 — H(l/m^). 

More formally, if the algorithm proceeds to round t, then ||s|| > e'^. By the definition of set Gi, we have 
Si, < (1 + ^)s *2 < • • • < (1 + < 6 - 54 . Hence, > (sijef > ^{pt-i)/{me‘^) > (e'/(v^e))2, 

so the surpluses of goods in Gi are similar up to a factor of e and bounded from below. Also we have 

(S4 - > {sije{m + 1))^ > <^{pt_^)/{e'^{m + ifm). 

Since we rely on an approximate demand oracle, the surpluses of goods in Gi might not change in a 
monotone fashion when increasing their prices. Nevertheless, we can relate the surplus in the beginning 
and the end of a round as follows. For every i E Gi, by Claim the surplus from exact demands 
satisfies s' < x • s*. Thus 

s' < s' + /i < XSi + < x{Si + /i) + /U = XSi + (1 + x)/i. 

Since x < 1 + , it holds that (1 + x)/r < 3^ < Si/{R 2 m^). This means the increase within a round is 

bounded by s' < (1 + )si. Since we do not touch the price of any good j E G 2 , the WGS property 
implies for exact demands s'- > Sj. Hence s'- > Sj — 2/r. 

Now, in order to bound the change of ‘h(pt), we consider s' according to Gi and G 2 . We distinguish 
two cases. 


Case 1: max{s( | i E G 2 } < /r Intuitively, in this case the algorithm has decreased the surplus of some 
good in Gi to approximately 0 (recall that /r is sufficiently small). This decrease alone brings down the 
potential function by a factor of 1 — H(l/m). All other surpluses will cause an increase by a factor of 
at most 1 + 0(l/m^). 

More formally. Claim [A~4l gives us ^ < min{s' | i E Gi} < Ifi. Hence, the contribution of goods of 
Gi to <h(pt) can be upper bounded by 


k 





fc-i 


<1: 



— 1 

R2nrP ) 


2 


+ 4^. 


Furthermore, for every z E G 2 , if —< s' < we have sf < and if s( < —/r, by the WGS 

property of the market, we know Sj < s' < s' + /i < 0. Thus, since s'- > Sj — 2/r, 


E ' 

3^02 


^ E 

3 &G 2 

s'.<-inP 


{Sj 2/i) < 


E 

teGa 
' <-Tra^ 


1 + 


R2nrP 


Hence, the contribution of goods of G 2 to <h(pt) can be upper bounded by 

-/ 1 ^ 

1 + 


E 

j=k+l 


< max ■ 

j=k+l 


R2m'- 


sl , 


< i: (' 

j^k+1 


R2m^ 


si 
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Combining the two parts 


$(pO = ^st < E(i + 


i=i 






§1 + (m® + 49 )/x2 


^ ^ R2m^ 


< 

< ii- 


R^rn^ 

1 

2e^m 


($(pt_i) - SiJ + (m^ + 49)/x 


^’(Pt-i) 


1 - - 
{ 

$(Pi-i) 


196 


R\w7 R\m}^ 


where the last inequality holds for any m > 2 with sufficiently large constants i?i, i ?2 • 


Case 2: max{s^ I * £ G 2 } > ^ Intuitively, in this case the gap between surpluses in Gi and G 2 decreases 
to approximately 0. Below we show that the closing this gap yields a decrease of the potential function 
<I> by a factor of 1 — n(l/m^). All other surpluses will increase by a factor of at most 1 + 0(l/m^). In 
combination, it turns out that <I> will decrease by a factor of 1 — n(I/m^). 

More formally, in this case min{s'j | i G Gi} > max{s'- | j G G 2 }. Let = min{s'j | i G Gi} and 

SG 2 = max{s'- I j G G 2 }. For every i G Gi, let s'^ = x'si — Si where x' = (1 + ), and for every j G G 2 , 

let s'j = Sj — 2fj, + 5j. Hence Si,Sj > 0 for all i,j- Further, we have | J2ih and | < mfi , 

hence 


and 


Now we have 


E = E - E 

i^Gi 

> E + E = E + E ~ 

ieGi jeG2 ieGi jeG2 

> ^ Sj — 4m/U 
j&G2 

1 1 

iSGi 


^(pt) = E = E Gsi - Sif + {sj - 2 ^+ SjY^ 

i i^G-i j^G2 


E + E - 2m)' 

yiGGi ieG2 

+ E “ 2Ai) 

jeG2 


Y -2^i + Sj)- Y^ Si{x'si - <5i) - E 

\j^G2 i^G\ j i^G\ 


(7) 


< X 


^4>(pt_i) - 4^ ^ Sj +4mfj,‘^ + SG 2 Y ~ E E 


ieGa 


iGG2 *eGi 


i^G\ 


+ (Sifc+l 2^) 'Y^ Sj 


j^G2 


< x''^^{Tpt-i) - + SsGim^i +- {Si^ - Y, Si 


( 8 ) 


(9) 


i^G\ 


< x'^4>(pt_i) + 4 to^^ + 24m^/r - -(si^ - 


( 10 ) 
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^ \ R 2 'm? R^rrfi Riui'^ 4e^(m + l)^m/ P‘ ^ 

Here [7] can be derived by expanding the quadratic formula and appropriately reorganizing the terms. 
For the step from ([T]) to ([5]), in the first bracket we overestimate the quadratic terms of s into x'^ 
and 1^21 by m. In the second bracket, we return to s'^ and §'■, which in turn are bounded correctly using 
sgi for all i G Gi and sg 2 for all j G G 2 . For the final two terms in ([7]) and ([5]) we use the definition of 
and and the fact that x' > 1. For the step from ([5D to Q, for the first bracket of ([5]) we use 
Sj > Sm for every j E G 2 . For the second bracket of ([5]) we note sci > SG 2 the difference between 
the sums of J-terms is bounded by 4m^ as noted above. By the same argument, we can transform the 
last two terms of ([5]) as shown. Note that we simply drop —2fj.J2j£G2 ^ From Q to (fTUll we use 

the fact that every surplus is bounded by 2m in its absolute value by Claim ETTl For the last term we 
use the bound for J2ieGi noted above. From (fTOll to (fTT|) . we just replace pL by its definition and use 
the bound <h(pt_i) > and (s^ — > ‘h(pt_i)/(e^(m + l)^m) as shown above. 

Finally, using sufficiently large constants Ri,R 2 , the multiplicative term in (jllj] can be decreased to 
strictly less than 1 for every m > 2. The final expression in (|12l) captures the asymptotics and proves the 
lemma. □ 

The previous lemma shows a decrease in the potential only for rounds in which the x determined by 
binary search is rather small. Lemma 12.21 continues to hold here and bounds the number of rounds with 
a larger value of x. The following variant differs only in the constant R 2 , and its proof is literally the 
same as for Lemma 12.21 

Lemma A.2. During a run of Alg-WGS, there ean be only 0{m‘^Di) many rounds that end with 

a: > 1 + p ^ 3 . 

Finally, we can assemble the properties to show the main result. 

Theorem 12.11 For any market that satisfies Assumptions \KJ\ and and for any e > 0, Alg-WGS 
returns the priee veetor of an {1 + e)-approximate market equilibrium in time polynomial in the input size 
and log(l/e). 


( 11 ) 

( 12 ) 


Proof. Let xt be the value of x we find in round t of Alg-WGS. First because at least one price will 
increase bv a factor of xt in round t. bv Claim IA.3I we have Ut Xt < 2™-^!. At the end of round t, if 
Xt > ^ + fot s = max{si I Si G 0{pt-i)} and s' = max{sj | s* E 0{pt)}, then by Claim ?? 

we have ‘h(pt) < ms'^ < mxls^ < mx‘l^{pt-i). Moreover, by Lemma IA.2I there will be at most 
0{m^Di) such rounds. Hence the total increase of <h(pt) in these rounds will be no more than a factor 

of nx.>i+i/iJ2m3 mx? < 

For all other rounds, we have x < 1 -|- , and by Lemma lA.il the potential function is decreased 

by a factor of 1/(1 — Therefore the total number of rounds before <h(pt) < will be at most 


/ / 2 1 \ 

O I logi/(i_n(^))-- 1 = O ( DimGogm-tmMog-j 


In each round, the number of queries to the oracle is no more than 0(log A) = 0{Di + D 2 +\ogm + \ogj). 
We conclude that the total number of queries during the algorithm is 0{(^Dimf \ogm + mA log {Di + 
D 2 -|- logm -|- log f)), which is a polynomial in the input size and log(l/e). □ 
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B Exact Equilibrium for Exchange Markets with Spending Constraint 
Utilities 

Here we show how to obtain an exact market equilibrium in exchange spending constraint markets. Us¬ 
ing the Alg-Spending-Exact framework we convert the approximate market equilibrium obtained in 
Section 13.21 into an exact equilibrium. To achieve this, we rely on the full information of the spending 
constraint utilities. This step is an extension of the technique developed in [M] for the linear exchange 
markets. However there are several challenges due to the much more involved setting of spending con¬ 
straint utilities, where the allocated partitions make the remaining budgets of agents and the values of 
goods dependent on too many parameters. In the following we present how to handle the extra complex¬ 
ity of the problem and this result resolves an open question of [24j of finding an exact polynomial time 
algorithm for exchange markets with spending constraint utilities. 

Let p be the price vector of an (1 -|- e)-approximate equilibrium. We first construct a bipartite graph 
EG'(p) = (A U G,E) where the edges E is a union of equality edges in EG{p), edges due to positive 
endowments, and edges due to allocated segments. The main idea here is to construct a set of components 
of agents and goods such that there is no interaction across the components. In addition, we want at 
least one good with price 1 in every such component. To achieve this latter condition, whenever there is 
a component C of EG'{p) without a good with price 1, we raise the prices of goods in C by a common 
factor X > 1 until a new equality edge appears. By Assumption 13.11 a new equality edge will always 
emerge in during this procedure because prices of goods in G are increasing, which makes goods outside 
G more and more attractive to the agents in G. 

The next lemma shows that the updated price vector after the price increase still remains a (1 -t- e)- 
approximate equilibrium. 

Lemma B.l. The price vector p at the end of while loop in Alg-Spending-Exagt is a (1 -|- e )- 
approximate equilibrium. 

Proof. Note that we increase prices of goods in a component G when each good has price greater than 
one. Corollary 13.11 implies that the surplus of each good in G is zero. Hence, the old allocation will still be 
feasible after the price change and the surpluses remain the same. Therefore the updated p will remain 
a (1 -|- e)-approximate equilibrium. □ 

At this stage we can assume that each component of EG'{p) has a good with price 1. We then work 
on each component of EG'{p) separately. We assume for convenience that EG'{p) is a single component. 

Next we set up a system of linear equations in price variables of the form Ap = b, and show that 
the matrix A has full rank. Finally, we will show that by perturbing the vector b slightly we can get 
an exact equilibrium. Consider the components of EG{p), i.e., after removing edges due to endowment 
and allocated segments from EG'{p). Let Gi,... ,Gk be the set of components of EG{p). In each 
Gi,l < I < K, all goods are connected with each other through a set of equality edges. Whenever there 
are two current segments (i, j, k) and (i,/, k') of the same agent i, we have the following relation between 
the prices of goods j and j': 


^ijkPj' ^ij'k'Pj- (1^) 

This implies that for each component G/, |G/ C G| — 1 of these equations are linearly independent, and 
there is essentially one free price variable. Further, since there is no money flow across components with 
respect to the current allocations, we have the budget balance condition for each component; 

Remaining worth of goods - remaining budgets of agents (after allocated segments) 

= sum of surpluses 
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Algorithm 5. Alg-Spending-Exact 

Input : Exchange market with a set A of agents and a set G of goods; Wij,Uijk, Bijk are market 
parameters as defined in Section 11.11 
Output: Prices p of an exact market equilibrium 


1 m •(— |G|; n •(— |A|; L •(— total bit length of all input parameters; e 

2 p ■(— (1 + e)-approximate equilibrium using Alg-Spending with Alg-Spending-Rounding 

3 s •(— 0(p). If s = (0,0,..., 0) then return p. 

4 EG{p) •(— (undirected) equality graph at prices p /*as defined at the beginning of Section [ST/ 

5 E ■(— {{i,j,k) I {i,j,k) is an allocated segment} 

6 EG'{p) <r- EG{p) U {{i,j) I Wij > 0} U {{i,j) I B for some k}. 

7 while EG'{p) contains a connected component C that does not has a good with price 1 do 


8 

9 

10 


Find the smallest x > 1 such that EG{p) C EG(Update(p, x, C)). 
p ^ Update(p, X, G). 

Recompute EG{p) and EG'{p). 


11 /* Wlog we assume EG'{p) consists of only one connected component. If there are more than one, 
then apply the procedure below individually to each component */ 

12 Let Cl,..., Gk be the connected components of EG{p) 

13 Set up the following system of linear equations in price variable 


1. p' = 1 for a good i whose price is 1 


2. For each component Ci,l < I < K 

(a) |C;| — 1 linearly independent equations of the form Uijkp'j/ = Uij/k'p'j, where (i, j, fc) and 
(i,/, k') are the current segments. 

(b) J2jeCinGP'j - J2j RijP'j = 0, where 

Rlj = J^ieCiCiA “ J2{i,j',k)eF-,j'^CinG Rij'k) + '^ij ,k)&F-,j'&Gir\G Rij'k 

14 p' G- the solution of above system of equations 

15 returu p' 


For component C;, the condition reads 


(Pj ^ Bijk^^Wijipji) ^ C^'^ij'Pj' ^ Bijk'^^Wijipji) — ^ £j, 

i&CiDG {i,j,k)£F j' ieGitlA j' (i,j,k)£F j' j&GinG 

where E is the set of allocated segments. Rearranging the above equation, we get: 
j&GiHG j&G ieC,nG 

where Rij = 51 Wij{l - 51 Bipk) + Y Y 

ieCiOA {i,j',k)eF-,j'^GinG i^GiDA {i,j',k)eF-,j'eGinG 

Each Rij is a rational number with denominator at most 2^^, where L is the total bit length of all input 
parameters Wij,Uijk and B^jk- 

Lemma B.2. For every 1 < I < K and j S C/ H C, 0 < Rij < 1. For every j G G, J2i Bij = 1- 

Proof. For each good j, we have 'Zi Wij = 1. Further both E(i,i',fc)eF;i'^c,nG Bij'k and 'E{i,j',k}eF-,j'eGi Rij'k 
take values in [0,1], hence the hrst claim of the lemma follows. For the second claim, Yli Bij = ^ij — 
1 . □ 
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Let M be the coefficient matrix of the system of equations (|14p . Then 


Mij = 


i-Rij ifjsQnG 
—Rij, otherwise 


From Lemma [B.21 each column j of M has at most one positive entry, namely Mij where j G C; H G, and 
each column of M sums to zero. There are in total L equations of type (1141) . one for each component. 
Next we eliminate the equations of type (I13p . Then there will be only one price variable per component. 
We designate a representative good for each component, say good I for Ci. Then each price pj in Ci is a 
constant multiple of price pi of good 1. Let pj = ajpi, where aj is a rational number whose numerator 
and denominator are products of at most m Uijk’s. Now we can rewrite the budget balance equation (fT^ 
for Cl in terms of L price variables as follows: 


Pi ^j-Y.Pi' Y. Y 

jeC,nG v j&Ci,r\G j&Glf^G 


Let Ti = Y.j&CinG^j^^ii' = EjeC,,nG ei = Y.j&Gir\G^i- The above equation becomes: 


PiTi - '^puSiv = £i 

v 


(15) 


Let N be the coefficient matrix of the system of equations (I15p . Then 


Nw 


Ti - Su if 1 = 1' 
—Sill, otherwise 


Since both Ti and Sw are rational numbers with denominator at most 2™^ and 2^*”^ respectively, each 
Niii is a rational number with denominator at most 2^™^. 

Lemma B.3. 0 < Nu for every 1; Nir < 0 for every I V; and J2i ^iv = 0 fof every I'. 

Proof. The proof essentially follows using Lemma IB. 2 1 The first two claims are straightforward, and for 
the last claim we have 

Y^ii' = Y Y = Y (^jY^ij= Y = 

i i jeC,/nG jeG;/nG i jeC,/nG 


□ 

Since there is a good i with price 1, we assume with loss of generality that good i belongs to component 
K, hence cxkPk = 1- The next lemma is an adaptation of a result of |24j . 

Lemma B.4. The K equations eonsisting of equation (USD for components 1,..., K — 1 and the equation 
OiKPK = 1 are linearly independent. 

Proof. Let N' be the coefficient matrix of this system of equations. It is easy to check that N' is the 
same as N except that = 0, l<i<iL — 1. Assume by contradiction that there is a non-zero vector 
a = (ai,..., Ok) such that a'^N' = 0. Let a/g be the entry in {oi,..., ax-i} that has the largest absolute 
value, and with loss of generality we assume that ai^ > 0 and the first K' entries of a are equal to ai^, 
i.e., ai = ■ ■ ■ = ax' = azo- 
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For each I < K' we have 

0 = ^ ahNhi + ax -0 

l<h<K 

= “Zo X! ^hi-aif^NKi+ X! {o,h - aio)Nfii 
l<h<K K'<h<K 

= -cliqNki + ^ {ah - aiQ)Nhi 

K'<h<K 

Using Lemma lB.3[ the above implies that Nhi = 0 for K' < h < K and I < K'. Next we show that 
Nih = 0 for 1 <l < K' and K' < h < K as well. By summing up the equations (I15|) for 1 < / < K', we 
get 

^ ez = ^ Niipi + Y^ NihPh 

1<K' 1<K' h l<K'-,l^h 

= Y + X! X! ^IhPh + X! X! ^IhPh 

1<K' h<K'l<K';l^h h>K'1<K' 

= Yj ^hhPh + X! X! ^IhPh + X! X! ^IhPh 

h<K' h<K'l^h h>K'l<K' 

= Y. PhY.^lh+ Y Y. NlhPh 

h<K> I h>K'1<K' 

= Y, Y, ^IhPh- 

h>K' 1<K' 

Since Ph > 1 for every h, if some Nih is non-zero for I < K' and h > K', then the right-hand side is at 
most —1/2^™-^, which is a contradiction. It implies that Nih 7 ^ 0 iff both I and h are either less than 
or equal to K' or larger than K'. Let Ai and Gi denote the set of agents and goods of components 
Cl,..., Ck’ respectively. Let ^2 = ^ \ and G 2 = G \ Gi. We can further conclude that 

• Wij = 0 for every i ^ Ai, j E G 2 , otherwise Nih / 0 for Z < iL' and h > K'. Similarly Wij = 0 for 
every i E ^ 2 , 3 E Gi. 

• Agents in Ai have no allocated goods in G 2 , otherwise budget balance equations of components 
Cl, I > K' will have some non-zero PjA ^3 ^ and that will make Nih 7 ^ 0 for I > K' and 
h < K'. Similarly agents in A 2 have no allocated goods in Gi. 

This is impossible since we have assumed that EG'{p) consists of one single component. Therefore N' 
must have full rank. □ 

Overall, we have established that equations of (fT^ . (fT^ and Pi = I are linearly independent. We can 
write this system in the matrix form as Ap = b, where A is invertible and all entries are rational numbers 
with common denominator at most 2 ^^. 

Consider the system Ap' = b' for a price vector p', where b' is a unit vector with a one in the row 
corresponding to the equation pi = 1. Next we show that p' gives an exact equilibrium. For that we need 
to show the following: 

(a) Equality edges with respect to p' and p are the same. This will imply that all allocated segments 
remain allocated. 

(b) (s, A U G U t) is a min-cut in N{p'). Combining the equation(2b) in step [12] of Alg-Spending- 
Exact, this will imply that there is a feasible allocation on current segments which gives surplus 
of each good. 
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By Cramer’s rule and Lemma lB.2[ the solution of Ap' = 6' is a vector of rational numbers with common 
denominator D < That is, all p'j^ are of form qi/D, where qi,D are integers. Since 

11^1 —1^11 < 2e, we have \p'i—pi\ < 2eD for every i. Let e' = 2D'^e, then \Dpi — qi\ = D\pi—p'^ < 2eD'^ = e'. 
For part (a), suppose UijkPj' < Uij/k'Pj then 

‘Uijkqj' ^ Uijki^Dpj' T £ ) ^ DUij'k'Pj T Uijk£ ^ '^ij'k'qj T i^ij'k' T Uijk')£ 

^ '^ij'k^qj T f 

Since both Uijkqj' and Uij/k'qj are integers, we have Uijkp'j/ < Uij/k'p’j- This implies that all equality edges 
with respect to p' and p are the same. That further implies that all allocated segments remain allocated. 
For the part (b), consider the network A^(p^) with respect to prices ph Recall that in iV(p') 

• The capacity of edge from source node s to agent i is J2j' Wij'p'j,{l — J2j,k-.{i,j,k}eF ^ijk)- 

• The capacity of a MBB edge from agent i to good j is J2j' Wij'p'y. 

Since pCs are all rational numbers with a common denominator D < all capacities in 

A^(p') are rational numbers with a common denominator no more than D^. Also because \p[ —pi\< 2eZ), 
the capacity of each edge e in A^(p) is at most the capacity of e in N{p') plus AeD. 

Let c be the capacity of cut (s, A U G U t) in A'(p')- Suppose there is a min cut in A^(p') with value 
less than c. Then that value is at most c — l/Z?^. This same cut in A^(p) will have value at most 
c — I/D'^ + (m + n + mn)4:Ds. Also the capacity of the cut (s, A U G U t) in A'(p) is at least c — AneD. 
Therefore the total surplus of goods in A^(p) is at least 

c — AneD — (c —— (m + n + mn)A£D = — 4(m + 2n + mn)£D > e, 

V D‘^ J 

which is a contradiction. Hence condition (b) also holds. 

We conclude with our main theorem. 

Theorem 13.11 For any spending constraint exchange market satisfying Assumption \S. 1[ Alg-Spending- 
Exact returns the price vector of a market equilibrium in time polynomial in m and L. 

Proof. By Lemma [3. 6 1 we know that we arrive at a (1 + e)-approximate equilibrium in time polynomial in 
m, L, and log(l/e). For the exact equilibrium, Alg-Spending-Exagt solves a system of iL < m linear 
equations whose entries are polynomially bounded in m and L, hence can be done in polynomial time. 
Further, log(l/e) is a polynomial of m and L as £ = l/m^™'2^'” Since M is a polynomial of m, L, 
binary search and Rounding run in polynomial time in each round of the framework. Hence, the total 
running time of the algorithm is polynomial in m and L. □ 
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